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Abstract 

The topological susceptibility of the vacuum in quantum chromodynamics has been simulated 

numerically using the Asqtad improved staggered fermion formalism. At nonzero lattice spacing the 

residual fermion doublers (fermion "tastes" ) in the staggered fermion formalism give contributions 

to the susceptibility that deviate from conventional continuum chiral perturbation theory. In 

this brief report we estimate the taste-breaking artifact and compare it with results of recent 

J^j ■ simulations, finding that it accounts for roughly half of the scaling violation. 

CO 

^ ■ PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.39.Fe 
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I. INTRODUCTION 

The staggered fermion formalism reduces the number of fermion doublers to four species 
("tastes") per quark flavor. At zero quark mass in the continuum limit the taste symmetry 
group becomes SU(4)£,xSU(4) R xU(l)y, but at nonzero lattice spacing this symmetry is bro- 
ken by lattice- artifact, taste-changing interactions. Taste breaking effects modify predictions 
of chiral perturbation theory, complicating efforts to extrapolate lattice results simultane- 
ously to zero quark mass and zero lattice spacing. However, a recently developed modifi- 
cation to chiral perturbation theory, called "staggered chiral perturbation theory (S^PT)", 
permits a quantitative treatment of taste-breaking effects [1, 2]. Combining these modifica- 
tions with an improved staggered fermion formulation that substantially reduces the strength 
of taste breaking has proven highly effective in carrying out high-precision extrapolations of 
results from lattice simulations [3]. 

In this brief report we use staggered chiral perturbation theory to predict the topological 
susceptibility. The derivation is given in Sec. II and in Sec. Ill the resulting formula is 
compared with results of recent lattice simulations with improved staggered fermions. 

II. TOPOLOGICAL SUSCEPTIBILITY WITH TASTE SYMMETRY BREAKING 

We consider the low energy effective chiral Lagrangian for N flavors of staggered fermions, 
each with four tastes. The Lagrangian is written in terms of the (AN) 2 pseudoscalar meson 
fields (ftia, encapsulated in a unitary 4A x 4A matrix 

U = exp(iT ia (f) ia /f). (1) 

We use / ~ 130 MeV. We label the generators of the Lie algebra of U(4N) with a flavor 
index % and taste index a and write them as 

T ia = Ait a , (2) 

where the t a , the 16 taste matrices (Dirac gamma matrices plus the identity), are enumerated 

as 

6, £a*> i£iis = «£/*&, i^u = i£ M &/, $i = l (3) 

and the Aj are the N 2 generators of U(N). The t a are orthonormal: 

Tr(t a t b ) = 45 ab , (4) 



but we normalize the Aj so that 

Tr(A,A,) = 5 l3 . (5) 

Moreover, for convenience we choose a quark flavor basis in which the diagonal generators 
are enumerated first and are explicitly 

(Ai) jfc = SijSik for i = 1, . . . , N. (6) 

Also, for convenience, we define the flavor-singlet generator and fields 

A = / (7) 

N 

^oa = $>«,/>/#. (8) 

The Euclidean Lee-Sharpe S%PT Lagrangian [1], extended by Aubin and Bernard to N 
flavors [2], and including an explicit mass term for the anomalous flavor singlet field is 

C = ^Tr (d^d.u) - ^Tr[M(tf + U)] + ^ + j v{u) (g) 

where a is the lattice spacing, the trace is over all 4N taste and flavor indices, Ai is the 
taste-degenerate quark mass matrix 

Miajb = M i:j 5 ab = SijSabrrii, (10) 

and the taste-breaking potential, 

-v(u)=Y,QO u (ii) 

is a linear combination of the operators 

1 = Tr (T 0i5 [/T 0i5 [/t) (12) 

2V = ~ [Tr(T ^U)Tr(T 0t ^U) + h.c] (13) 

2A = ~ [Tr(To, At5 f/)Tr(To,5^) + h.c] (14) 

3 = ^[Tr(T ^UT 0;tl U)+h.c.] (15) 

4 = ~ [Tr(T 0iliB UT ofili U) + h.c] (16) 
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Tr(T 0i/ ^)Tr(To^t)j (17) 



o 5V = - 

5A = l - [Tr(r , M5 t/)Tr(r 0)5M f/t)] (18) 

6 = J2 T* (T Qtlu ,UT , vli rf) . (19) 



The partition function is, as usual 

Z= f[dU]exp(-f CdV). (20) 

The chiral condensate is defined by 

E 4 = lim 1^ (21) 

and is, as usual, 

S 4 = S = fif/2 (22) 

The bare meson masses are obtained by expanding the Lagrangian to quadratic order in 
the meson fields: 

£ = \ E(<^a) 2 + \ E AiajbMjb (23) 

ia iajb 

where the flavor and taste mixing matrix is 

A iaJb = A^dab + A\fS Ia 5 Ib + a^Babdij. (24) 

The continuum taste-singlet flavor-mixing matrix is block diagonal with an invariant sub- 
space % = 1, 2, . . . , N. This is the only subspace that contributes to the tree-level topological 
susceptibility. For i,j in this set we have 

Af = 2/iTr(MA i A i ) = 2 f xS ij m i (25) 

Af) = ml/N. (26) 

The flavor-singlet taste-mixing matrix B a b is defined by the small field expansion 

W) -\Y. BabMu- (27) 

i,a,b 

Without the anomalous mass term tuq and without taste breaking the bare flavor-neutral 
meson masses are given by the diagonal values in A^ '. 

m 2 ia = 2fmii for i = 1, 2, . . . , N. (28) 

To lowest order in the taste-breaking potential the pseudoscalar mass splittings were calcu- 
lated by Lee and Sharpefl]. With the Aubin and Bernard generalization [2] they are 

Am 2 ia = « 2 ^a, (29) 



where A a = B aa and 



A 5 == (30) 
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A M = --(67! + (73 + 3(74 + 3(76) (31) 



16 
A M5 = jj (Ci + 3C 3 + C7 4 + 3C 6 ) (32) 

16 
A^ = — (2(7 3 + 2C 4 + 4C 6 ) (33) 

Aj = ^(4C 3 + 4C7 4 ) . (34) 

To formulate the topological susceptibility in the presence of taste breaking, we recall 
the derivation, starting from the quark and gluon Lagrangian Cqcd- The partition function 
receives contributions from topological charge sectors v as follows [4]: 

Z = Y,Z» > ( 35 ) 

V 

where the microcanonical partition function is 

Z v = J dqdqdAd9exp(-i9v)exp ( I\Cq C d + id F F / Un 2 ]dv\ . (36) 

The theta term in the Lagrangian can be removed by a flavor (and taste) singlet {7a (1) chi- 
ral rotation, leaving a Lagrangian £qcd{@) with rotated quark fields. The lattice staggered 
fermion action breaks this symmetry through the quark mass term and through taste sym- 
metry breaking at nonzero lattice spacing a. To order a 2 in S^PT the symmetry breaking 
is described by the £/a(1) non-invariant taste-breaking term a 2 V(U). The C/a(1) rotation 
on the quark fields is modeled in the effective theory by the corresponding rotation on the 
meson matrix U: 

C[9) = ^Tr (djfldfjj) - ^Tr{M(e~ w / AN U^ + e ie ' m U)} + ^<f> 2 01 + a 2 V{e ie ' m U). (37) 
8 ^ ' 4 2 



A field redefinition <\>n — > <\>u — f9 /AN for i = 1, . . . N, implying O / —> 4>oi — f9/AvN, gives 
an alternate expression 

C(9) = ^Tr (d^dfjf) - !±£.Tr[M(rf + U)] + ^(0 O / - f9/AVN) 2 + a 2 V{U). (38) 

We work in the mean-field e regime in which the fields are space and time independent 
and take the limit mLV — > oo, and ttIqV — > oo i.e. asymptotically large mean square 



topological charge (v 2 ). In this limit we may use a saddle-point approximation to compute 
the microcanonical partition function. The saddle point occurs close to the identity in U and 
origin in 9, where we replace the Haar measure dll with the (AN) 2 dimensional Cartesian 
measure d<f> and expand the integrand to quadratic order in the meson field <fi and 9: 



Z u oc / dcf) d9 exp(—i9i/) exp 



V 



F((j)) + m' (<i>oi-fO/WN) 



(39) 



where F(4>) is the quadratic form 



J2 [4i ) *«6 + 2 Sa6«5i. 



>ia9jb- 



(40) 



ia,jb 



The only meson fields contributing to the susceptibility are the diagonal taste-singlet mesons 
4>u for % = 1,2, ... ,N. The corresponding subspace is invariant under the quadratic form. 
Thus the remaining fields can be integrated out immediately. The integration over theta 
gives 

r V 



Z u oc / ]J <i0j/exp(— i4VN(p 0I i//f) exp 

J i=l 



-F(^) - 16Nu 2 /(2Vfm 2 ) 



(41) 



and finally the integration over taste-singlet meson fields gives 



N 



Z u oc exp 
where the bare taste-shifted flavor-neutral masses are 



16v 2 / ^VfY.yK) - lQNu 2 /(2Vfm 2 ) 



(42) 



m 2 j = 2/imi + a 2 Aj. 



(43) 



The resulting topological susceptibility is 



X = (v 2 )/V 



f 2 /16 



Ef=i V^ 2 / + N/ml 



(44) 



When the N flavors can be grouped into p degenerate flavor multiplets with degeneracy n^ 
so that Y?k=i n k — N, the susceptibility is 

/ 2 /16 



X=(S)/V 



(45) 



E P k= in k (l/m 2 kI + l/m 2 ) 

Counting the taste symmetry as well, the continuum limit degeneracy of quark species 
is then An k for each distinct flavor. In numerical simulations of staggered fermions it is 
necessary to take fractional powers of the quark determinant to simulate a theory with the 
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physical count of 2 + 1 flavors, i.e., only two species of light u and d quarks of the same 
mass and a single s quark. The S^PT treats only the case in which quark multiplets come 
in multiples of four. Absent a rigorous treatment of fractional powers in S;\PT, we adopt 
the replica trick: assume that an analytic continuation in the flavor number n^ reproduces 
the effect of taking a fractional power of the determinant. That is, we replace 4n u d by 2 and 
An s by 1 in the expression above. The susceptibility for 2 + 1 flavors is then 

X ~ 1 + mljpmlj + 3m^/2mg (46) 

Of course, this result is valid only in leading order chiral perturbation theory, but at large 
quark mass it interpolates smoothly [5] to the quenched result predicted for pure gauge 
theory [6]. 

Xq = fm 2 /12. (47) 

At small m u d the susceptibility is dominated by the pion mass. Without taste breaking we 
recover the conventional formula 

X = f 2 m 2 j8. (48) 

We see that the effect of taste breaking at leading order is simply to replace the masses of the 
Goldstone mesons in the conventional expression with the masses of the corresponding taste- 
singlet non-Goldstone states. The squared masses of these states vary linearly with quark 
mass and have a nonzero intercept. Therefore, at nonzero lattice spacing, the topological 
susceptibility does not vanish at zero quark mass. However, in the continuum limit, the 
expected zero is recovered. 

III. COMPARISON WITH SIMULATION RESULTS 

The topological susceptibility has been determined recently for gauge configurations gen- 
erated with an improved staggered fermion action[7]. It was found that at nonzero lattice 
spacing the lattice susceptibility is higher than that predicted by chiral perturbation theory, 
but (at least within the available statistical errors and an extrapolation based on limited 
data) it is consistent with the continuum prediction. Two principal causes for disagreement 
at nonzero lattice spacing were proposed: (1) Dislocations in the gluon field could imitate 
small instantons and participate in screening of the topological charge. These and small 
instantons are erased in the smoothing (cooling) process, possibly leaving behind larger, 



unscreened instantons. The latter would then increase the smoothed susceptibility. (2) The 
would-be zero modes are not at zero because of taste-symmetry breaking, so are not fully 
effective in suppressing the topological charge at small quark mass. The latter effect is taken 
into account in Eq (46). 

In light of our result [Eq (46] we replot reanalyzed data of Ref. [7] as a function of the 
taste-singlet pion mass m^ jTq in Fig. 1(b). For comparison the same data is plotted as a 
function of the Goldstone pion mass rn^r^ in Fig. 1(a). The reanalysis is described in detail 
elsewhere [8]. In this work we are concerned with the taste-breaking artifact. 
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FIG. 1: Topological susceptibility vs pion mass squared in units of tq with improved (Asqtad) 
staggered quarks. Data is from Ref. [7], reanalyzed with a new variance reduction method and new 
0.09 fm data point [8]: (a) plotted as a function of the Goldstone pion mass squared (b) replotted 
as a function of the taste-singlet pion mass squared. The bursts show the continuum extrapolation 
at fixed abscissa in each case. The solid line shows the prediction of leading order continuum 
chiral perturbation theory [Eq (46)] without the U(l) mass term ntQ. The dashed line includes ttiq, 
determined through Eq (47) from the simulated quenched susceptibility [8]. 

We see that correcting for staggered fermion artifacts removes roughly half the scaling 
violation in this range of pion masses. The rest can be attributed to the definition of 
the topological charge density operator, including smoothing [7]. It is plausible that those 
corrections decrease as 0{a 2 ). We assume this form in carrying out the rough continuum 
extrapolation. The extrapolation to the continuum starts by fitting the measured points for 
each lattice spacing to the curve 1/x = cjm^ 1 + b and then performing an extrapolation 
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linearly in a 2 of the thus smoothed values of \. In the left panel the extrapolation was done 
at constant Goldstone pion mass m 2 ; in the right panel, at fixed m 2 j. Because staggered 
fermion artifacts contribute to scaling violations at fixed Goldstone pion mass, and those 
artifacts decrease as 0(a 2 a 2 ), rather than 0(a 2 ), the extrapolation in Fig 1(a) is somewhat 
less plausible than the extrapolation in Fig 1(b). 

To make further progress we need simulations at a smaller lattice spacing. 
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